The enigma of the 1^ = quantum Hall effect in graphene 
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We apply Laughlin's gauge argument to analyze the v = Q quantum Hall effect observed in 
graphene when the Fermi energy lies near the Dirac point, and conclude that this necessarily leads 
to divergent bulk longitudinal resistivity in the zero temperature thermodynamic limit. We further 
predict that in a Corbino geometry measurement, where edge transport and other mesoscopic effects 
are unimportant, one should find the longitudinal conductivity vanishing in all graphene samples 
which have an underlying v — quantized Hall effect. We argue that this u — graphene quantum 
Hall state is qualitatively similar to the high field insulating phase (also known as the Hall insulator) 
in the lowest Landau level of ordinary semiconductor two-dimensional electron systems. We establish 
the necessity of having a high magnetic field and high mobility samples for the observation of the 
divergent resistivity as arising from the existence of disorder-induced density inhomogeneity at the 
graphene Dirac point. 



A single two-dimensional (2D) layer of carbon atoms 
forming a honeycomb lattice, i.e., a graphene layer, has 
unusual physical properties attracting a great deal of 
current interest.1 Among its intriguing properties, the 
effective low-energy dispersion is linear in 2D momen- 
tum: E — ?iw|k|, where v, the graphene Fermi velocity, 
is a constant {v ~ c/300 where c is speed of light), and 
electron wave functions formally obey a Dirac-like con- 
tinuum equation with zero Dirac mass, rather than the 
Schrodinger's equation. Associated with this Dirac na- 
ture is the fact that the single-particle spectrum has a 
two-fold pseudo-spin (or valley) degeneracy, in addition 
to the usual (double) spin degeneracy. Thus the low- 
energy spectrum of graphene is made of particle or hole 
excitations near a double-cone Fermi surface; the apex of 
the cones are the Dirac points where electron and hole 
dispersions cross each other, or where the valence and 
conduction bands become degenerate. The combined 4- 
fold spin/pseudospin degeneracy gives rise to an emer- 
gent SU(4) symmetry, which is useful in analyzing the 
low-energy properties of graphene and plays a central role 
in our discussion below. Obviously this is a very unusual 
band structure, which has attracted much attention both 
theoretically and experimentally. 

The observation of quantum Hall effect (QHE)^'9, when 
an external, perpendicular magnetic field (B) is applied, 
provides the most compelling evidence for the 2D mass- 
less Dirac nature of electrons in graphene. In particular, 
such a system is predicted^i^ to support integer QHE 
with quantized values of Hall conductance given by: 



±gsgvin+l/2)e^/h, 



(1) 



with n = 0, 1, 2, • • • is an integer, and gs = gv — are re- 
spectively the spin and pseudospin/ valley degeneracies; 
the latter is inherent in the chiral, massless Dirac equa- 
tion describing 2D graphene. The half integer form in Eq. 



(HI) arises from the Berry phase associated with the pseu- 
dospin index, and the experimental observation of the se- 
quence predicted in the form of Eq. ([T]) is a direct reflec- 
tion of the massless chiral Dirac nature of the low-energy 
electronic states in graphene. Aside from the shift 1/2, 
Eq. ([T]) implies two additional peculiar characteristics of 
graphene QHE: (i) The ± sign in front of the quantized 
Hall conductance indicates the smooth transition from 
hole-like to electron-like carriers as Fermi energy moves 
through the Dirac point, as expected from the gapless 
Dirac spectrum; (ii) the degeneracy factor g = gsgv = 4 
is associated with the SU(4) symmetry discussed earlier; 
in the specific case of massless Dirac electrons moving in 
an (orbital) magnetic field, it is equivalent to the state- 
ment that each Landau level (LL) has a, g = gsgv = 4 
fold degeneracy, on top of the usual orbital degeneracy. 
The graphene LL energy: 



En = ±y/2nnv^\eB\ 



(2) 



with n = 0, 1,2,---, is yet another key difference be- 
tween Dirac spectrum as compared to ordinary electron 
LL spectrum: En = {n + \/2)fiujc with ujc = eB/{mc), 
where m is effective mass of the electron (or hole) in a 
given band. The spectrum of Eq. ^ has been observed 
directly in STM measurements.^"* We note that the strict 
2D nature of graphene (i.e., absence of any higher sub- 
bands due to the finite width of quantum wells in usual 
semiconductor materials) and the large inter LL energy 
separation for large B and small n has made it possi- 
ble to observe a room temperature quantized Hall resis- 
tance plateau in graphene; in particular, the axy = 2e^ /h 
plateau is observed with a quantization accuracy of 0.2% 
at T = aOOif and B = 45T.9 

To facilitate comparison with QHE in regular semicon- 
ductor systems, we rewrite Eq. ([T]) as 

a-^y = ve^ /h, (3) 
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with V = ±g{n + 1/2) and g = Qsgv = 4. For SU(4) 
symmetric ground states, we expect QHE at 



±4n + 2 = ±2,±6,±10, ■ 



(4) 



which are indeed the observed QHE sequence at modest 
magnetic fields; in fact, this sequence has been observed 
up to = 26 at B = 9T and T = 1.6K; this observation 
decisively reflects the chiral massless Dirac dispersion of 
graphene with a pseudospin Berry phase. 

The question we are interested in pursuing is whether 
the SU(4) symmetry can be broken in the ground state, 
leading to the splitting of the 4- fold degeneracy. This lift- 
ing has, in fact, been observed at higher magnetic field 
B > 2QT, where new QHE has been seen for ly — 0,±1, 
and ±4iiS Theoretically this issue has been discussed at 
great length in the literature fii^i^ Various sources of the 
splittings, including Coulomb interaction, electron spin 
Zeeman splitting, lattice effects not captured by Dirac 
equation, and electron-phonon coupling, have been iden- 
tified, but no consensus has been reached thus far. In- 
stead of adding to the discussion of the microscopic origin 
of the splitting, we take a phenomenological approach in 
this paper, and discuss the consequences of the splitting 
when present. 

To understand these anomalous quantum Hall (QH) 
states, we introduce spin and valley splittings to the 
Dirac LL spectrum: 
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= ±v/2n?i«2|es|/c± A»(B) ± A;'(B), (5) 



n = LLs 



n = 1 LLs 



where A", A" are the spin and valley splittings of the nth 
orbital LL respectively, which depend on both the LL in- 
dex and magnetic field B. The spectrum in the presence 
of such splittings is illustrated in Fig. 1. For definitive- 
ness we assume A^, < A^ <^ l-En+i ~ En\, which are rea- 
sonable assumptions for the LLs of interest in graphene. 

We first ask, assuming finite A^ and A^, on phenomeno- 
logical grounds, what the expected QHE sequence would 
be for the LL stricture of Eq. ^ and Fig. 1. Making the 
standard QHE assumption of the LL tails being occupied 
by localized states and the LL centers being extended, we 
conclude that, if all the 4- fold degenerate LLs are indeed 
split into distinct spin and valley split individual levels, 
we will have the sequence of QHE with 

i. = 0,±l,±2,... , (6) 

instead of the sequence of = zLAn+2 = ±2, ±6, ±10, • • • 
with SU(4) symmetry present in QH ground states. If we 
assume As > while A^ — 0, so the symmetry reduces 
from SU(4) to SU(2), the corresponding sequence would 
be 

i^ = 0,±2,±4, ••• . (7) 

The same would be true in the opposite case with A^ = 
while Ai, > 0. 



FIG. 1; Upper panel: Schematic single electron density of 
states in graphene when the SU(4) symmetry is completely 
broken. As and A„ are the spin and valley splittings respec- 
tively. The spin and valley resolved Landau levels (LLs) are 
further broadened by disorder, with states extended at the 
center of each LL and localized away from the centers. The 
dashed lines indicate the locations of the Fermi energies at 
the center of representative quantum Hall plateaus. Lower 
panel: Schematic single electron density of states in Si (100) 
inversion layer when the SU(4) symmetry is completely bro- 
ken. 

The quantization of axy is associated with the vanish- 
ing of longitudinal conductivity axx in QHE, as can eas- 
ily be seen from Laughlin's gauge arguments (and aug- 
mented by Halpcrin).^'^ Thus corresponding to any of the 
QH sequence discussed here, we also have 

axx = (8) 

on the QH plateaus, at T = 0. We emphasize that the 
quantization of (Jxy and vanishing of axx in the limit 
T — > are intrinsically tied through gauge invariance, 
as is most clearly demonstrated by the Laughlin gauge 
arguments; these properties in fact define QHE. 

The main difference between 2D graphene and ordi- 
nary 2D semiconductor systems with parabolic or non- 
relativistic electron dispersion (e.g.. Si MOSFETs, GaAs 
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heterostructures and quantum wells) is not the 4-folcl 
spin-pseudospin degeneracy, but the presence of the n = 
LL dX E = Q which has equal weight in the electron and 
hole bands. It is instructive to compare the LL structure 
of graphene with that of the Si(lOO) inversion layer^i^ 
where the LL energies are 

£:„ = (n + l/2){heB/mc), (9) 

with n = 0,l,2,---, while the QH sequence is given by 

zy = gn = 4,8, 12,--- . (10) 

We note that Si(lOO) MOSFET electrons also have an 
SU(4) spin-valley symmetry, and consequently, each LL 
is f/-fold degenerate with g = gsgv = 4; yet the SU(4) 
symmetric QH sequence above is different from that of 
graphene, Eq. ([4]). Here the spin degeneracy is the same 
as in graphene, but the valley degeneracy, gy = 2, arises 
from the bulk Si band structure which has six equivalent 
ellipsoid minima (i.e., "valleys") in the conduction band 
close to the edges of the Brillouin zone; this is in contrast 
to the graphene case, where the valley degeneracy is an 
inescapable consequence of its Dirac nature. In a strong 
magnetic field the spin and valley degeneracies of Si(lOO) 
MOSFET are also lifted, giving rise to the LL spectrum 
(c.f. Eq.® for graphene, and see Fig. 1): 

En = (n + l/2){heB/mc) ± Al'{B) ± A[^{B). (11) 

This would lead to the familiar integer QH sequence i> — 
1,2,3, Comparing with that of graphene Eq. ([6|) 
when SU(4) symmetry is completely broken, other than 
the negative v states corresponding to hole QH states, 
the biggest difference is the presence of i/ = QH state 
in graphene, which is the focus of this paper. 

It is instructive to discuss further the phenomenology 
of QHE in Si inversion layers in the context of under- 
standing graphene QHE since both systems in their pris- 
tine states have spin and pseudospin (i.e. valley) symme- 
tries. The most extensively studied Si QHE is in Si (100) 
inversion layers where the original QHE discovery was 
made by von Klitzing.J>S This 2D system, as mentioned 
above, has an SU(4) symmetry, which is found to be lifted 
in the low-lying LLs leading to = 1, 2, 3, 4, • • • quantiza- 
tion in the observed QHE p^y = h/{ve'^). We emphasize 
that in Si (100) 2D system, the valley splitting is much 
less than the spin splitting, and in fact, in higher or- 
bital LLs the valley splitting is typically not manifested in 
the observed QHE which follow the SU(2) sequence with 
V = 2,4, 6,8, • • • for larger orbital levels. This is again 
similar to graphene where valley splitting effects have so 
far seemed to have been observed only in the lowest LL. 
Similar to graphene, the origin of valley splitting in Si 
(100) system is still not theoretically well- understood, ^'^ 
and there is not much direct experimental evidence for 
valley splitting in the absence of a magnetic field. It 
is conceivable that in graphene the valley splitting is 
strongly enhanced by interaction effects in a magnetic 
field leading to the full lifting of the SU(4) symmetry. 



The phenomenological parallel between graphene and Si 
QHE with respect to the ground state SU(4) symmetry 
breaking points toward the possible importance of inter- 
action effects playing a dominant role in lifting the valley 
degeneracies in both cases. Finally, we point out that in 
Si (111) inversion layers, the ground state valley degen- 
eracy is six, leading to an effective SU(12) ground sate 
symmetry. Random uniaxial interfacial stress can lift this 
SU(12) symmetry converting it to an SU(4) ground state 
symmetry by splitting the six-fold valley degeneracy into 
two lower valleys and four upper valleys. Therefore, Si 
(111) system becomes equivalent to the Si (100) SU(4) 
situation, with the possible further lifting of the SU(4) 
symmetry in high fields at the lowest LLs. QHE has been 
seen in Si (111) experimentally.-^ 

In Si MOSFET and other non-relativistic 2D electron 
systems, for sufficiently high B (where all electrons are 
in the n — Q LL with filling factor v <^ 1), the system 
enters an insulating regime with longitudinal resistivity 
Pxx — > oo as T — > 0, while the Hall resistivity pxy tak- 
ing its classical (and roughly temperature-independent) 
value. It is worth noting that experimentally, one nor- 
mally uses the Hall bar geometry to directly measure 
Pxx and Pxy (actually to be more precise, longitudinal 
and Hall resistances Rxx and Rxy, and convert them to 
Pxx and Pxy with appropriate geometrical factors); axx 
and then in turn be obtained through the tensor 

inversion relations: 

(Txx = PxxUpIx + Ply)\ <^xy = Pxy/iplx + Ply) ■ (12) 

It is clear that in this high B insulating phase (often 
called the Hall insulator) one has axy and axx 
in the limit T ^ 0, for a non-zero range of B\ this, by 
definition, corresponds to a = QH state! But in 
that context this high B insulating phase is not referred 
to as a QH state, mainly because no plateau is seen in 
Pxy, which is what is measured directly, despite the fact 
that there is a plateau in axy, as we argued above. This 
relates to our discussion about the important difference 
between QH states with v = Q and v ^ Q. 
The tensor inversions of Eqs. are 

Pxx ^ Crxx/{(^lx+ (^lyY Pxy ^ Crxy/{(ylx+ (^ly)- (13) 

As a result for any QH state with 7^ 0, we have 

Pxx = Q; Pxy = l/a-xy = h/{iye'^). (14) 

These equations are deceptively simple, and form the ba- 
sis of studying QHE - the deceptive dichotomy here, of- 
ten lost in the literature, is that theoretically the fun- 
damental quantities are axxi= 0) and axy{= ve^jh), 
whereas experimentally the fundamental (or more specif- 
ically, the measured) quantities are pxx and pxy 

The a, p dichotomy has not been an issue in QHE stud- 
ied in semiconductor based 2D systems because there one 
always has 

Oxy ^ Oxx^^ (15) 
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in the quantization regime, leading to Eqs. (|14p : thus 
the longitudinal transport coefRcients vanish and Hall 
coefficients quantized, in both a and p. 

In graphene, however, the n — LL resides precisely at 
the particle- hole symmetric point of E — 0, which allows 
for a (Txy — QHE when the SU(4) symmetry is broken, 
as indeed observed experimentally i"'^°'^^i^°'^^i^^'^'^ Putting 
a^y = in Eqs. IT^ we get 



P: 



xy 



0; 



(16) 



But, the quantization of Hall conductivity, even when the 
quantized value is zero, demands axx — *■ as T — > 0, and 
consequently we have for the v = graphene QHE 



Px 



(17) 



As discussed earlier, this is very similar to what hap- 
pens in the high field insulating phase (or Hall insulator 
phase) in ordinary 2D electron gas with parabolic band 
structure. The parallel between them is actually quite 
complete. In particular, the experimental Hall insulator 
phase has pxx going to infinity and Pxy given by the clas- 
sical formula, pxy = B/Nec, where N is the 2D carrier 
density in the system. As we discussed above, an inver- 
sion of the p tensor leads to the Hall insulator conductiv- 
ity (Jxx = and axy = 0, which is formally the same as 
the v — QH conductivity values around the Dirac point. 
Now we note that the corresponding graphene Hall insu- 
lator phase must necessarily have pxy = since the 2D 
carrier density at the Dirac point (i.e. at the zero energy 
charge neutrality point), separating the electron and hole 
bands, is by definition zero (i.e. iV = in the classical 
Hall resistivity formula). This shows that the emergent 
divergent longitudinal resistivity at the graphene Dirac 
point could equally well be considered a putative Hall 
insulator phase, which in this context, is equivalent to 
the zero Fermi energy v ~ Q QHE in graphene with 
o'xx = and pxy ~ 0. The only (qualitative) differ- 
ence between graphene and 2D semiconductor systems is 
that the graphene Dirac point is known to be dominated 
by density inhomogeneities associated with the electron- 
hole puddleS j^^i^^ i^'^^'^^ which lead to considerable den- 
sity fluctuations around the average (expected) zero den- 
sity at the Dirac point. This means that there will be 
considerable fluctuations around the expected Pxy — 
value in the graphene Hall insulator phase. If these fluc- 
tuations are large, then neither the divergent resistivity 
nor the v — Q graphene QH (or the equivalent Hall in- 
sulator phase) would be observable, indicating the need 
for very high mobility samples where these fluctuations 
are suppressed^^ as well as very high magnetic fields so 
that the magnetic length is well below the typical puddle 
size. We emphasize that in this bulk picture that we are 
proposing there is no finite-T phase transition in the sys- 
tem, only a possible T = quantum phase transition to 
the V = Q phase. At finite T and in the presence of strong 
edge effects (i.e. mesoscopic samples), there may be sub- 
stantial modification of this underlying picture, but we 





v=0 








(b) 



^ r 



FIG. 2: (a) Schematic illustration of formation of domains 
with different f's due to electron and hole puddles, and cur- 
rent carrying states flowing along the domain walls, (b) A 
building block (saddle point) of an appropriate network model 
describing these domain wall states. 



think that it is unlikely that there will be any universal 
physics in such a situation and details (such as the mech- 
anism producing the SU(4) symmetry breaking) will play 
crucial roles. 

Based on very general Laughlin-type gauge considera- 
tions applied to bulk 2D graphene QHE, we thus arrive 
at the result summarized below: If the SU(4) spin/valley 
symmetry of graphene ground state is lifted, then there 
is a. V — Q QH state with axy = when the Fermi en- 
ergy is at (and around) the charge-neutral Dirac point, 
with a divergent bulk longitudinal resistivity pxx oo 
when T 0; this is a necessary consequence of hav- 
ing a = QH plateau. We thus believe the re- 
cent observation22i2ii22i2^ of divergent resistivity at the 
graphene neutral point in a high magnetic field is a di- 
rect manifestation of the SU(4) symmetry- broken v — 
QHE. 

In the following we make a few more comments about 
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this conclusion, in which some specific issues related to 
experimental studies will be addressed. 

(i) All our arguments presented above are about the 
zero-temperature bulk physics in the thermodynamic 
limit; technically, this implies that the thermodynamic 
limit is taken before any other limits, in particular, the 
zero temperature limit T — s- 0. As a result meso- 
scopic effects like edge transport play no role in our 
considerations. It is conceivable however, that small 
graphene samples used in experiments may be in the 
mesoscopic regime at sufficiently low temperatures, and 
edges may play an important or even a dominant role 
in transport; transport theories based on edge chan- 
nels have been developed, "^^'-^ It is also conceivable 
that difference o-^°'^^i^°'^^i^^'^'^i'^° observed in the transport 
properties near the Dirac point may also be due to the 
differences in edge properties of different samples. This 
leads us to a straightforward prediction: If one measures 
the bulk conductivity directly using the Corbino geome- 
try, one should get axx ^ in all samples exhibiting the 
z/ = QHE, due to the absence of edge contribution to 
transport. 

(ii) Our bulk arguments are completely independent of 
the details of how the SU(4) symmetry is broken. It could 
be due to a spontaneous symmetry breaking driven by 
Coulomb interaction, or an explicit symmetry breaking 
induced by Zeeman splitting or intervalley scattering in- 
duced valley splitting. On the other hand the mesoscopic 
effects discussed above are sensitive to the specific man- 
ner in which the symmetry is broken, see, e.g., Ref. Isol 
and on whether the valley-splitting or the spin-splitting 
is larger in magnitude. 

(iii) Our bulk consideration suggests that pxx ~* oo 
does not reflect a bulk phase transition driven by B at 
finite T, as it is a property of the underlying QH phase. 
There can be a quantum phase transition at T = 0, 
driven by B, into the 1^ = QH phase*^ Temperature 
can, however, be very important in the role played by the 
above mentioned mesoscopic effects. 

(iv) It is known,— as discussed above, that due to 
the presence of disorder arising from charged impuri- 
ties and the lack of screening due to vanishing density 
of states at the Dirac point, neutral graphene samples 
are quite inhomogeneous, and electron and hole puddles 
form at i? = 0. Due to such density inhomogeneity, it 
is quite possible that the observation of the 1^ = QHE 
in graphene necessitates not only the complete lifting of 
the SU(4) symmetry in the n — graphene LL, but 
also a sufficiently high magnetic field. This is because 
in the presence of B, these electron and hole puddles 
will turn into domains of different (integer) z^'s, with cur- 
rent carrying states propagating along the domain walls; 
see Fig. 2a. Such domains and domain walls on top 
of the 1/ = background (or percolating domain), when 
present (and especially if large), can dominate &mZA; trans- 
port, thus suppressing the divergent pxx- The situation 



is very different for v 0, where the background or per- 
colating domain has a non-zero Hall conductance, thus 
dominating bulk transport. As B increases, we expect 
the symmetry breaking splittings and/or A„ (that 
give rise to the v = Q QHE) to increase correspondingly 
This will lead to suppression and eventual elimination 
of these domains. We thus conjecture the critical field 
observed in Refs. ,20,. 21. beyond which very strong in- 
sulating behavior kicks in corresponds to the field that 
leads to such suppression. This is clearly consistent with 
the observation that the more disordered the sample (as 
measured by mobility), the higher the critical field , as 
the electron hole puddles or domains are due to disor- 
der. Theoretically, one can formulate a network model 
to describe transport through these domain wall states; 
see Fig 2b. It turns out the appropriate network model 
takes the same form as that of a system with random but 
zero average magnetic flux.'^^ This system was argued to 
support a Kosterliz-Thouless (KT) type metal-insulator 
transitioni^ A number of numerical works indeed find ev- 
idence of such a transition;^ although it was argued later 
on— that actually all states are localized, albeit with 
extremely large localization length at weak disorder, giv- 
ing rise to an apparent transition in finite size numerical 
studies. Such an apparent transition may be related to 
the apparent KT-like transition seen in Refs. [20ll2ll . We 
therefore suggest that scanning probe measurements22i2S 
be carried out in graphene samples under QHE condi- 
tions to directly probe the electron-hole puddles and the 
associated domains which may be dominating the exper- 
imental QH physics around the Dirac point. 

We conclude by asserting that theoretical considera- 
tions of the V — QHE in graphene using Laughlin's 
gauge argument lead to a divergent bulk pxx associated 
with the quantization ofaxy, when the Fermi energy is at 
the Dirac point. This is an intrinsic bulk property of the 
system, which can be viewed both as a quantized Hall liq- 
uid and a Hall insulator. In our consideration edges play 
no role, and there does not need to be a bulk phase tran- 
sition associated with the the divergence of pxx ■ A clear 
and simple prediction is that all systems that exhibit a 
quantization of axy = will have axx — > in a Corbino 
disk geometry measurement. We also suggest that the 
experimentally observed magnetic field-driven insulating 
behavior is due to the suppression of domains or elec- 
tron and hole puddles by field. These domains should be 
directly observable in scanning probe measurements. 
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ures. This work is supported by Microsoft Station Q 
(SDS), DARPA-QuEST (SDS), NSF-NRI (SDS), and 
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out while the authors were visiting Kavli Institute for 
Theoretical Physics (KITP). The work at KITP is sup- 
ported in part by National Science Foundation grant No. 
PHY-0551164. 
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